ANNALES 
UNIVERSITATIS MARIAE CURIE-SKLODOWSKA 
LUBLIN — POLONIA 
VOL. XV, 10 SECTIO AA 1960 


Z Katedry Fizyki Teoretycznej Wydziału Mat.-Fiz -Chem. UMCS 
Kierownik: doc. dr Włodzimierz Urbański 


Maksymilian PIŁAT 


Drgania wielowarstwowej ciekłej kuli z rdzeniem 


Ko.sreóaana MHOTOCIOŃHOTO ZKMĄKOTO Wapa € TBEPĄBIM CEPĄEHHUKOM 


The Oscillations of the Many-Layered Liquid Sphere with a Core 


$ 1. Teorię małych drgań grawitacyjnych kropli nieściśliwej cieczy ideal- 
nej podał Lord Kelvin [1]. Lamb uogólnił te wyniki na przypadek 
jednej warstwy cieczy pokrywającej sztywny rdzeń sferyczny [1]. Lord 
Rayleigh- opracował teorię drgań kapilarnych kropli takiej cieczy. 
W teorii jądra atomowego u podstaw modelu kroplowego [3], a także 
uogólnionego tzw. kolektywnego, rozpatruje się drgania kropli cieczy 
idealnej i nieściśliwej o stałym rozkładzie gęstości masy i gęstości ładunku 
tj. p= const., q = const., gdzie p oznacza gęstość masy, q — gęstość obje- 
tościową ładunku elektrycznego. Uwzględnia się przy tym napięcie po- 
wierzchniowe kropli obliczone z energii wiązania jądra. 

Celem mojej pracy jest uogólnienie wspomnianych wyżej wyników 
na następujące przypadki: 

1) drgań n warstw różnych nieściśliwych cieczy idealnych, naładowa- 
nych elektrycznie z różnymi gęstościami objętościowymi, koncentrycznie 
pokrywających sztywny rdzeń sferyczny, 

2) drgań kropli cieczy o ciągłym rozkładzie gęstości masy i gęstości 
ładunku. 

W przypadku pierwszym częstości własne drgań dane są wzorem: 


gdzie: 
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Ze wzorów tych otrzymuję jako przypadki szczególne wyniki podane 
przez Lorda Kelvina i Lamba [1] oraz wzór na częstości własne 
drgań kropli pokrytej warstwą innej cieczy. 

W przypadku drugim wzór na częstości własne otrzymuję z wyżej 
wypisanych wzorów przez przejście graniczne, gdy liczba warstw rośnie 
nieograniczenie. 

$ 2. Rozpatrzmy drgania swobodne n warstw cieczy idealnych i nie- 
ściśliwych, koncentrycznie pokrywających sztywny rdzeń sferyczny. Niech 
R, Rp ..., R, oznaczają promienie powierzchni tych warstw przed od- 
kształceniem w kolejności licząc od warstwy zewnętrznej do rdzenia. Po 
odkształceniu powierzchnie warstw dane bębą równaniami r; =T„ (©, 9). 
Przyrosty Ar,=r, R,, rozwijam w szereg funkcji kulistych Y,, (©, ç): 


TE RZ [1 Te: ai Yu (9, p] (2, 1) 
gdzie ay) = ayn (t) sa współczynnikami zależnymi od czasu, a sumowanie 
rozciąga się po r od—h do FÀ i po h od 0 do oo, k=0,1, ..., n. Oczy- 


wiście r„=R„ i 4” =0. 
Zakładamy, że ruch cieczy jest bezwirowy, co łącznie z warunkiem 
nieściśliwości daje równanie dla potencjału prędkości ®,: 


Vq D Er bp 2 11 0. 


Drgania wielowarstwowej ciekłej kuli z rdzeniem 
Rozwiązaniem tego równania jest 


Bw 
p = > (a> r. + =) Y,, (8, 9) dla R„<r<R;y_ 
Me 


Potencjał ©, musi spełniać następujące warunki: 


W1) ©, jest ciągły w cieczy. 


p 
W2) v,= LA jest ciagla w cieczy. 
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Z warunku W3° wynika: 


A pati 
By, = Ay. Tq Rn 


Z W3 mamy dla k=0: 
A, RE" — A+ 1)B,, = RP ao 


p 
Łącząc wzory (2,5) i (2,6) otrzymuję: 
R 43 
Aj = DEET > DFI ao 
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Z W3 i wzory (2,5) wynika: 
1-3 
Rk * (k) 


Ay = X(R2T Rey an dla k=1, ...,n 
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(2, 2) 


(2, 3) 


(2, 4) 


(2, 5) 


(2, 6) 


(2, 7) 


(2, 8) 


(2, 9) 
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następnie ze wzorów (2,7) oraz (2,9) otrzymuję 
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Po scałkowaniu wzoru (2,10) otrzymuję 


Re = Re R 143 
(bt) —- Se (o) (k) 
"u Rat RH Ra). te (2, 11) 


sa stałymi określającymi początkowe odkształcenia powierzchni 
warstw. W rozpatrywanym przeze mnie przypadku ce =Q 


§ 3. Energie kinetyczną obliczam ze wzoru 


gdzie c” 
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Wyrażając grad P we współrzędnych biegunowych otrzymuję wzór 
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obliczam z (3,1) i (3,3): 


2 
„| 


Ry_, 
Py (24+ 1)/A,,| | 34041) (241) Bella dr = 
p X 1 1 1 
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k=1 

Po wstawieniu do powyższego wzoru 4,, i B,, ze wzorów (2,7) i (2,8) 

i po prostych przekształceniach otrzymuję: 
1 
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§ 4. Zmianę potencjału elektrostatycznego 59 związanego z odkształ- 
ceniem warstw naładowanej cieczy obliczam, traktując przemieszczenie 
ładunków w pobliżu powierzchni jako wytworzenie warstw powierzchnio- 
wych równych (q,,, q) Ar. 8% spełnia równanie Laplace'a, można je 
więc zapisać w postaci: 


żę, = ) (eer r+Q, up Yig (4,1) 


dp 
dla R,,<r<R,_,, gdzie k=0,1,..., n+1, R 
są współczynnikami zależnymi od czasu. 


n+ 
09 musi spełniać następujące warunki: 
W4) 09 jest skończony dla r=0 i lim òç = 0. 
T->oo 
W5) òq jest ciągły w całej przestrzeni. 


W6) Składowa normalna pola elektrycznego posiada skoki na gra- 
nicach warstw: 
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Z warunku W4 wynika, ze QU =Qrt?=o 
Z W5 wynika: 


Qa = Q*) R2H = CA BL = Qa” 
dy 


Warunek W6 daje: 


X (Q= QE) RO Q Ae 1) (Q?- ary =- Áx CR - da) Ry a (k) 


Łącząc wzory (4,3) i (4,4) otrzymuję: 


` 
— 47 (deti dy) Bia a 
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Zauważywszy, że (2a  =0 otrzymuję z (4,5): 


k 
QV _ 47 Gir % KO 
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© z (2,11) do (4,6) otrzymuję: 
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Wstawiajac œ 


gdzie K=O as NH. g. =U. 


Warto zauważyć, ze Qy""=Q. 


Ze wzorów (4,3) i (4,4) wynika, że: 
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Wstawiajac do (4, 9) ay. ze wzoru (2,11) otrzymuje: 
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§ 5. Zmianę elektrostatycznej energii kropli przy deformacji jej kształ- 
tu można przedstawić w następującej postaci: 


1 L ea 1 
AE =+ | aaeactt feta” att | 89-49 de = 


n+1 n—l = 
1 É 1 ; 
-+) Í 9,49, de ++ > Í e-dag de +4 > Ag,41°AG, dr (5,1) 
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gdzie Aq, = q, Ar, ò (r À R,), przy czym f Aq, dr=0 (5, 2) 


Korzystajac ze wzorów (2,1), (2,11), (4,1), (4,7), (4,10) i (5,2) otrzy- 
muje w rozpatrywanym przyblizeniu, ze: 
n=1 
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W celu obliczenia catki = á o. Aq* dt przedstawiam potencjał ?Q;., 


w cienkiej powierzchniowej warstwie (gdzie A9, 0) wzorem przybliżonym: 


bi 4/37 (R, = Rip) Vi b 4/3 Ta (R; — R.) 
r 2 a 
k+1 ` ` R. + Ar. Ry, de "kp 


RE 


przy czym R„,,=0. Z powyższego wzoru oraz ze wzorów (2,1), (21, 1) 
i (5,2) otrzymuję: 
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Zatem zmiana energii elektrostatycznej wyraża się wzorem: 


iag Ż, CH ao | (5, 6) 


n A k á 
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$ 6. Zmianę energii powierzchniowej przy odkształceniu kropli daje 
wzór [3]: 


NNEDYET 6 -00F2/a, F 
M 


(6, 1) 
gdzie 5 oznacza stałą napięcia powierzchniowego, a R promień kropli 
nie odkształconej. 


Dla n warstw analogiczny wzór będzie: 


n—1 
6. R Á 
= ` ` mrm a e 
k=o dye 


wstawiajac nastepnie ay ze wzoru (2,11) do (6,1) otrzymuję: 


n—1 : ` 
|: mai — U 2 R 2).46 
CY =A- NA+) > A mA pT) (R, (6, 3) 
R =R Ry 
k=0 n 
§ 7. Całkowita energia drgań wyrazi się wzorem: 


PAE 


Ga `) (T. 1) 
Ay 


gdzie C. =C! + c? (7, 2) 


(6, 1) 


gdzie: 
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Częstości własne drgań tej wielowarstwowej kropli dane są wzorem: 
„PS zb ve 3) 


§ 8. Przypadki szczególne: 

1. Kropla cieczy jednorodnej. Kładąc n=1. R, = R, = 0 otrzymamy ze 
wzorów przeze mnie wyprowadzonych wzory na częstości własne drgań 
swobodnych kropli cieczy stosowany w teorii jądra atomowego (3): 


87q° A — 1) _ PR, 


— 2 = ^ — Á 
C, = oR? (A — 1) @+2) 3R FI) 


2. Jedna warstwa cieczy na rdzeniu stałym. Kładąc we wzorach wy- 
prowadzonych w tej pracy n=1, R, =a, R, =R, =b, q =— Gp” i ozna- 
czając przez p, średnią gęstość rdzenia i warstwy cieczy oraz g = 4/3%G Pod, 
otrzymujemy po prostych przekształceniach: 


el 


P\g 
2 2 2+1 AET ) 
a+n(4) +(2] Fak 


a | 
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Wzór ten jest indentyczny ze wzorem podanym przez Lamba [1]. 

3. Kropla pokryta warstwą innej cieczy. Kładąc we wzorach (3,5), 
(6,3) i (5,7) n=2, R,=R,=0 otrzymamy wzór na częstości drgań kropli 
pokrytej jedną warstwą innej cieczy: 


C 
o =o gdzie C,=CY ke? 
a 


8m a> (A — 1) pH 4 ej 
"UA Ba 2 0 = 
C "aj 3 (2) + 1) (9; 91 % + 94) + 47a % alan 3 a 


3 
C? = (A—1) +42) a” [s+ (5) | 


5 22+1 
B, = le, — (P — P.) (2) | 


a ib są jak poprzednio promieniami warstw, tzn. R, =a, R, =b. 


$ 9. W celu obliczenia częstości drgań ciekłej kuli o promieniu a zciąg- 
lym rozkładem gęstości masy p= p (r) i gęstości ładunku 99 (r) zawie- 
rającej sztywny rdzeń sferyczny o promieniu b przechodzę we wzorach 
(3,5) i (5,7) z n do nieskończoności przy ustalonych R =a i R„=b. 
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Uwzględniam przy tym skoki gęstości na powierzchniach cieczy r=a 
oraz r=b. Otrzymane wzory mają postać następującą: 


= a 
URE 44 ala) _ @9 ) 8 a` 9 (a) 22+1 dQ 
ss < sy Ce: 13) tye |) dr «Aż 


' 22+1 
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a 
qs xa À qr” dr jest średnią gęstością ładunku w kuli. 
a 


o 
t 


W rozpatrywanym przypadku, napięcia powierzchniowe między war- 
stwami znikają w granicy i pozostaje tylko napięcie powierzchni zew- 
netrznej 5. Wtedy ze wzoru (6,1) otrzymuję: 


P= og Sf (1.1) A+2). 


Wzór na częstości drgań będzie taki, jak we wszystkich poprzednich 
przypadkach: 


= 
w? ar gdzie C,=C?+C?. 
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PESIOME 


PaccMOTpeHbl KOJIEÓaHMA N CJIOEB ZKMJĄKOCTK, IIOKPBIBAIOLNMX KOHIIEH- 
Tpuueckui TBepĄbii ceprqye4Huk. IlpeËxnojaraeM, YTO ZXUJĘKOCTE UĄCAJIBHA, 
HeCZKMMAEMA, HE MPOBOAMT SJIEKTPUHAECTBA, M YTO IIJIOTHOCTM Macc M 3a- 
PANOB B CJIOAX IIOCTOAHHBI. 

YureHo TAKKE IIOBEDXHOCTHOE HaTAZKeHME. BBIUMCJEHBI UACTOTBI COÓ- 
CTBEHHbIX KoJleóaHnii. Ilocpe4cTBoM riepexoqa K npełejry roJrydeHBI 
<bopMyJIbl uacror COÓCTBEHHBIX KoOJleÓaHui NIA CJIYUAA, KOTĄA IJIOTHO- 
CTM Macc M 3ApAĄOB ABJIAIOTCA HeEHpepbIBHbBIMM CbyHKUMAMM paccTOAHMA 
TOUKM OT WeHTpa mapa. 


SUMMARY 


The author investigated the oscillations of a n layer liquid, covering 
a rigid core concentrically. The liquid is assumed ideal. non-compressible 
and non-conducting, tne mass and electric charge densities being con- 
stant. The surface tensions have also been taken in account. 

The aim of the paper is to calculate the frequency of the self-oscil- 
lations. By means of a limiting process one obtains formulas for the self- 
oscillations in the case when the mass and electric charge densities are 
continuous functions of the distance of the point from the centre of the 
sphere. 
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Fliissiges NH,SCN*nNH; als Lösungsmittel. Teil I. Löslichkeit der 
Mettalle und einiger anorganischer Verbindungen. 


4 K. Zagórski: Hydrolityczne badania śluzu Inu. 
Hydrolytische Untersuchungen des Leinsammenschleimes, 


5 J. Czajka: Zmiana współczynników dó/de i dó/dT koloidalnego roztworu 
: białka w czasie termicznej denaturacji. 
Die. Veränderung der Koeffizienten dó/de und dó/dT der kolloidalen 
Eiweisslósung während der termischen Denaturation. 


6. D. Stachórska: Szybkość kondensacji pary przesyconej. II Kondensacja 
na jonach. 


The Rate of Condensation of Supersaturated Vapour.- II Condensation 
on Ions, 


7. J. Skierczyńska: Kilka uwag na temat pomiarów oporu właściwego Ge. 
Some Remarks on Resistivity Measurements of Germanium. 


8. J. Czajka: Studia nad wpływem stężenia, temperatury i czasu na napię- 
cie powierzchniowe koloidów hydrofilnych. 
Studien iiber den Einfluss der Konzentration, Temperatur und Zeit 
auf die Oberflichespannung der hydrophylen Kolloiden. 


9. S. Wieluński: Gerat zur Demonstration und Untersuchung der Zentri- 
fugal und Corioliskraft. 


Przyrząd do demonstracji i badań z siłami odśrodkową i Coriolisa. 
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